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Abstract. We study thermal entanglement in some low-dimensional Heisenberg 
models. It is found that in each model there is a critical temperature above which 
thermal entanglement is absent. 



1. Introduction 

Entanglement jlj plays an important role in quantum computation and quantum 
information processing. With appropriate coding, a system of interacting spins, such 
as described by a Heisenberg hamiltonian, can be used to model a solid-state quantum 
computer. It is therefore of some significance to study thermal entanglement in 
Heisenberg models. We find that for each model there is a corresponding critical 
temperature for transition to the entanglement regime, and the entanglement only 
occurs below this critical temperature. 

2. Measures of entanglement 

A pure state described by the wave function j^*) is non- entangled if it can be factorized 
as l^*) = l^*!) (8) 1^*2) • Otherwise, it is entangled. A typical example of an entangled 
state is the Bell state for a bipartite system of two qubits: 

^(101) -|10)) (1) 

For such a bipartite system the most popular entanglement measure is the 
entanglement of formation. For a pure state the entanglement of formation is defined 
as the reduced entropy of either subsystem|21. 

For the two-qubit system one can use concurrence 3 as a measure of the 
entanglement. Let pi2 be the density matrix of the pair which may represent either a 
pure or a mixed state. The concurrence corresponding to the density matrix is defined 
as 

C12 = max {Ai - A2 - A3 - A4, 0} , (2) 
where the quantities A^ are the square roots of the eigenvalues of the operator 

£'12 = (8cr|)pi2K '^cr|) (3) 
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in descending order. The eigenvalues of gi2 are real and non-negative even though gi2 
is not necessarily Hermitian. The entanglement of formation is a monotonic function 
of the concurrence, whose values range from zero, for an non-entangled state, to one, 
for a maximally entangled state. 

3. Heisenberg models 

The general A^-qubit Heisenberg XYZ model in a magnetic field B is described by the 
Hamiltonian 



n— 1 n— 1 

where we assume cyclic boundary conditions iV -|- 1 = 1 . The Gibbs state of a system 
in thermodynamic equilibrium is represented by the density operator 

p{T) = exp(-i//fcT)/Z, (5) 

where Z = tr[exp(— iZ/fcT)] is the partition function, k is Boltzmann's constant which 
we henceforth take equal to 1, and T is the temperature. 

As p{T) represents a thermal state, the entanglement in the state is called thermal 
entanglement. At T = 0, p(0) represents the ground state which is pure for the 
non-degenerate case and mixed for the degenerate case. The ground state may be 
entangled. At T = oo, p(oo) is a completely random mixture and cannot be entangled. 

4. Thermal entanglement in the 2-site Heisenberg model 

The density matrix can be obtained 4 as 

cosh( J/T) - sinh( J/T) 
- sinh( J/T) cosh( J/T) 

where A — (2 cosh( J/T) + 2 cosh(i?/T))~^, and the concurrence 

f sinh( J/T) - 1 \ 

^ = '""ncosh(J/T)+cosh(B/T)'°|- 
As the denominator is always positive, the entanglement condition is 

sinh( J/T) - 1 > or T < 1.134J. (8) 
from which we conclude that 

• There is a critical temperature Tc ~ 1.134 J. The thermal state is entangled when 
T<Te. 

• The critical temperature is independent of the magnetic field B. 

• Entanglement occurs only for the antiferromagnetic case (J > 0). 



p{T) = A 



(6) 
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(9) 



5. Thermal entanglement in the 3-site Heisenberg model 

We now consider pairwise entanglement in the 3-site Heisenberg model with uniform 
magnetic field, and also in the presence of a magnetic impurity The reduced density 
matrix of two sites can be written as 

f u \ 

w y 
y w 
V « J 

The concurrence may be readily obtained as 

C= ^max{|y|-V^,0}, (10) 

In the case of a uniform magnetic field, the entanglement between any two sites is 
the same due to cyclic symmetry. We therefore need only consider the entanglement 
between sites 1 and 2. Then 

u{B) = vi-B) = ^e'^^ + ie'3^(2z + z'^) 

w = cosh.{PB){2z + z-^) 
y = cosh(/3B)(z"2 - z) 

Z = 2cosh(3/3B) + 2cosh(/3B)(2z + z-2). (11) 

where (z = exp(/3J)). 

If _B = 0, one can easily find that the sites 1 and 2 are entangled if and only if 

2\z-^ ~ z\-3-2z- z-^ > (12) 

from which we conclude that 

• There is no entanglement when J > 0; 

• Entanglement occurs when J < and T < Tc, where the critical temperature is 
given by -1.27J = 1.27|J|. 

• The maximal concurrence is 1/3, which occurs for T — > 0. 

Fig.l plots the concurrence against r for different B. From these graphs we see 
that there exists a critical temperature above which the entanglement vanishes. It 
is also noteworthy that the critical temperature increases as the magnetic field B 
increases. 

We now consider the case of a single impurity field on the third site; thus 
Bi — B2 — and B3 — BJ > 0. In this case the cyclic symmetry is violated 
and we have to consider the entanglement between sites 1 and 2, and between sites 1 
and 3, separately. 

Fig. 2 plots the concurrence C12 and C13 against scaled temperature r — kt/\J\ 
for different magnetic fields B. From Fig. 2 (a) we see that when the magnetic field 
is located at the third site both the antiferromagnetic and ferromagnetic cases are 
entangled in the range < t < Tc, where the critical temperature Tc depends on 
B. Fig. 2(a) also suggests that the concurrence C12 tends to 1, namely that the 
(1,2) entanglement becomes maximal, when t for large enough i?, in both the 
antiferromagnetic and ferromagnetic cases. 

In contrast to the (1,2) case, the entanglement between sites 1 and 3 increases to 
a maximum with increasing B and then decreases. The lower the r, the smaller the 
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Figure 1. Concurrence as a function of T for different magnetic fields B = l(solid 
line), 3/2(daslied line), and 2 (circle point line). 




Figure 2. Concurrence C12 Cisagainst r for different B. For antiferromagnetic 
case (dotted line), B = 10. 
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B at which the concurrence reaches its maximum value. For smaller B, entanglement 
occurs only in the ferromagnetic case (J < 0), while for large enough B (e.g. i? = 10 in 
our units), weak entanglement occurs in both the antiferromagnetic and ferromagnetic 
cases. 

[1] Schrodinger E 1935 Naturwissenschaften 23 807 
[2] Werner R F 1989 Phys.Rev.A 40 4277 

[3] Hill S and Wootters W K 1997 Phys.Rev. Lett. 78 5022; Wootters W K 1998 Phys.Rev. Lett. 
80 2245; Coffman V, Kundu J and Wootters W K 2000 Phys.Rev. A 61 052306. 

[4] Wang X 2001 Phys.Rev. A 64 012313 Wang X 2001 Phys. Lett. A 281101 
[5] Wang X, Fu H and Solomon A I 2001 J. Phys. A 34 11307. 



